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We show that a simple model for the propagation of a 
rumor on a small-world network exhibits critical behavior at 
a finite randomness of the underlying graph. The transition 
occurs between a regime where the rumor "dies" in a small 
neighborhood of its origin, and a regime where it spreads over 
a finite fraction of the whole population. Critical exponents 
are evaluated through finite-size scaling analysis, and the de- 
pendence of the critical randomness with the network connec- 
tivity is studied. The behavior of this system as a function of 
the network randomness bears noticeable similarities with an 
epidemiological model reported recently [M. Kuperman and 
G. Abramson, Phys. Rev. Lett. 86, 2909 (2001)], in spite of 
substantial differences in the respective dynamical rules. 



Small-world networks have been introduced as an in- 
terpolation between ordered and random graphs, to cap- 
ture two specific features of real neural, social, and tech- 
nological networks On one hand, they have a rel- 
atively large clustering coefficient, i.e. the probability 
that two neighbors of a given vertex are in turn mu- 
tual neighbors is high, as in ordered graphs. On the 
other, the typical separation between two vertices is much 
smaller than the total number of vertices, as in random 
graphs. Whereas the geometrical properties of small- 
world networks have been studied in detail |p]-|l0|, less 
attention has been paid to the dynamical properties re- 
sulting from their partially random structure |p|,p^|^^^. 
As explained in more detail below, the structure of small- 
world networks is parametrized by the randomness p, 
where p = and p = \ correspond to fully ordered and 
random neworks, respectively. It has been shown that, 
in asymptotically large small- world networks, statistical 
geometrical properties — such as the mean distance be- 
tween vertices — exhibit qualitatively the same behavior 
for any p > 0. Here, we study an epidemic-like propaga- 
tion process on a small-world network which, in contrast 
with geometrical properties, shows a transition between 
two qualitatively different dynamical regimes at a finite 
value of p. Since this model addresses a social process 
and small-world networks are a presumably realistic rep- 
resentation of social networks, the analysis could be rel- 
evant to the description of threshold phenomena in real 
societies. 

Our model [ p^|Jl^ consists of an A^-element popula- 
tion where, at each time step, each element adopts one 
of three possible states. By analogy with epidemiological 



SIR models these states are called susceptible (S), 
infected (I), and refractory (R). The evolution proceeds 
as follows. At each time step a randomly chosen infected 
element i contacts another element j. Then, (i) if j is in 
the susceptible state, it becomes infected; (ii) if, on the 
contrary, j is infected or refractory, i becomes refractory. 
These rules are better interpreted in the frame of a rumor 
spreading process, where S-elements have not heard the 
rumor yet, I-elements have heard the rumor and are will- 
ing to transmit it, and R-elements have lost their interest 
in the rumor and do not transmit it. Initially, only one 
element is infected and the remaining N — 1 elements are 
susceptible. During the first stage of the evolution, the 
number of I-elements increases. Since this also implies 
a growth of the R-population, the contacts of I-elements 
between themselves and with R-elements become more 
frequent. After a while, in consequence, the I-population 
begins to decline. Eventually, it vanishes and the evo- 
lution stops. At the end, elements — now in the re- 
fractory state — have been infected at some stage during 
the evolution. Numerical simulations and analytical re- 
sults show that, generally, Nji < N. Therefore, there 
is a fraction of the population that never hears the ru- 
mor. In the original version of this model, each I-element 
is allowed to contact at random any other element of 
the population. It has been proven that, in such situ- 
ation, the ratio NfijN approaches a well-defined limit, 
Nji/N = 0.796 . . ., for asymptotically large values of A'' 

Here, in contrast, we assume that the elements are 
situated at the vertices of a small-world network, and 
contacts can be established between linked elements 
only. The small-world network is constructed from a 
one-dimensional ordered network with periodic bound- 
ary conditions — a ring — where each node is linked to its 
2K nearest neighbors, i.e. to the K nearest neighbors 
clockwise and counterclockwise [p],p|jl^]. Then, each of 
the K clockwise connections of each node i is rewired 
with probability p to a randomly chosen node j, not be- 
longing to the neighborhood of z. A short-cut between 
two otherwise distant regions is thus created. Double and 
multiple links are forbidden, and realizations where the 
small- world network becomes disconnected are discarded. 
As advanced above, the parameter p measures the ran- 
domness of the resulting small- world network. Note that, 
independently of the value of p, the average number of 
links per site is always 2K. We have performed series of 
lO'^ to 10^ numerical realizations of the model for several 
values of p, N, and K. At each realization, the small- 
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world network was generated anew and the evolution was 
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FIG. 1. Frequency distribution of the number Nr of 
R-elements at the end of the evolution, for K = 2 and (a) 
p = 0.05, (b) p = 0.19, and (c) p = 0.3. Different symbols 
correspond to N = 10^ (crosses), A'' — 10"^ (empty dots), and 
A'^ = 10^ (full dots). In (b), the dotted straight line has slope 
— 1.5. The dotted curve in (c) is a schematic representation 
of the A'^-independent profile observed for small Nr. Note 
carefully the different scales of the three plots. 

In the first place, we have studied the distribution of 
the number of R-elements at the end of the evolu- 
tion for K — 2. Figure |^ shows the normalized fre- 
quency f{Nfi) obtained from series of 10^ realizations 
for selected values of p and N. For small randomness 
{p — 0.05, Fig. |l]a) the distribution is approximately ex- 
ponential and does not depend on N. In this regime, 
the rumor "dies" after a few time steps in a small neigh- 
borhood of its origin, due to the lost of interest of the 
highly interconnected elements close to the initially in- 
fected site. Therefore, the size N of the whole population 
is irrelevant to the value of iV/j. The situation is consid- 
erably different for relatively large randomness [p = 0.3, 
Fig. I^c). Here the distribution /(A^_r) is bimodal, with a 
maximum close to Nfi = (not shown in the figure) and 
an additional bump for larger N^. Near Nn = 0, the 
frequency is independent of iV, as in the case of small 



randomness. Contributions to this zone of the distri- 
bution come from realizations where propagation ceases 
before a short-cut is reached. In contrast, the additional 
structure is strongly dependent on N. The position of its 
maximum, in fact, grows linearly, as N^'^^ « 0.25A^. In 
a typical realization contributing to this zone of the dis- 
tribution, many contacts occur through short-cuts and a 
finite portion of the population becomes infected. The 
intermediate regime, just before the large-iV/j structure 
begins to build up, is illustrated in Fig. |l|b for p = 0.19. 
The frequency follows here a power law, f{N]i) ^ N'^" 
with a w 1.5, over a substantial interval. This interval 
is limited by above by a smooth cut-off, which shifts to 
larger values of Nji as N^, with (3 ^ 0.5. 
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FIG. 2. Average fraction r — (Nr/N) of refractory ele- 
ments at the end of the evolution as a function of the ran- 
domness p on a small- world network with K — 2, and for 
N = 10^ (dots), N = 10* (dashed line), A^ = 10^ (dotted 
line), and A' = 10^ (full line). The insert shows a close-up in 
the transition zone. The arrow indicates the critical value of 
p, determined as explained in the text. 

The appearance of a well-defined power-law distribu- 
tion is a clue to the critical-phenomenon nature of the 
transition from the regime where the rumor remains lo- 
calized to the regime where it spreads over a finite frac- 
tion of the population. To characterize the transition, we 
choose as an order parameter the average fraction Nn/N 
of R-elements at the end of the process, 

N 

r = {NnjN) = ^ Nnf{NR), (1) 
NR=a 

which can be straightforwardly calculated from the nu- 
merical results. In the small-randomness regime, where 
f{Nji) is independent of N, we expect r ^ N~^, so that 
r —^ as N —^ oo. For large p, in contrast, the pres- 
ence of the large- A^fl^ maximum in /(A^i?) should provide 
a finite contribution to r even for asymptotically large 
populations. These features are verified in numerical re- 
alizations. Figure ^ shows r as a function of p, calculated 
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in series of 10^ to 10^ realizations for N ranging from 10'^ 
to 10^ and K — 2. A well-defined transition at a finite 
randomness pc ~ 0.2 is apparent for large N. The insert 
of Fig. H shows a close-up of the main plot near pc- The 
smoothness of the curves, observed even for the largest 
values of N, suggest that the critical exponent 7 associ- 
ated with r ~ |p — Pcl''^ just above the transition, should 
be larger than unity. 

To determine the values of Pc and 7 we apply finite- 
size scaling analysis Q . Results are presented in Fig. ||. 
The insert shows a plot of the fraction r as a function of 
N for several values of the randomness p. With p — 0.17 
and 0.18 we find the subcritical behavior expected for 
p < Pc, i.e. r ^ N^'^. For p = 0.20 and 0.21, on the 
other hand, r saturates to a finite value for very large N. 
Between these two regimes, withll p = 0.19, r decreases 
as r - N-P, with p = 0.78 ± 0.02. We identify this 
intermediate regime with the critical point, and thus get 
Pc = 0.19 ±0.01. 
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FIG. 3. Data collapse in the finite-size scaling of the frac- 
tion r as a function of p—pc, Eq. (^, with u = 0.36. Diflterent 
symbols correspond to different values of A*', ranging from 10'' 
to lO'^. The insert shows r as a function of A*' for five values of 
the randomness: from below, p = 0.17, 0.18, 0.19, 0.20, and 
0.21. The dashed and the full line have slopes —1 and —0.78, 
respectively. 

Once Pc has been determined, the critical exponent 7 
is obtained from the finite-size scaling Ansatz p8| 



rN^-" ^F[{p-pc)N-'] 



(2) 



For p = Pc we have r — F{0)N^"''^ , so that = p = 
0.78 ±0.02. The best data collapse near the critical point 
is obtained for i/ — 0.36 ± 0.01, as shown in Fig. |3[ This 
yields 7 = 2.2 ±0.1. 

We have examined this model for other values of K, 
up to K = 10, and found the same kind of transition 
in all cases. The average fraction r as a function of p 
is shown in Figure ^ for = 10^ and several values of 
K, calculated over 10^ realizations. It is seen that the 
critical randomness Pc decreases as K grows. This is due 



to the increment in the number of long-range contacts per 
element, as a consequence of the higher connectivity of 
each site. On the other hand, the value of r at p = 1, ri, 
grows with K and approaches the level expected for the 
original version of the model, r* = 0.796 . . . (l|0. The 
insert in Fig. || displays the dependence of pc and of the 
difference r* — ri with K . Though this plot covers less 
that one order of magnitude in the if-axis, the results 
suggest power-law decays for both quantities. 
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FIG. 4. Average fraction r of refractory elements at the end 
of the evolution as a function of the randomness p for several 
values of K, and A'^ = 10^. The insert shows the critical 
randomness pc (full dots) and the difference r* — r\ (empty 
dots; see text for definitions of r* and ri) as a function of K. 

In summary, our numerical analysis shows that the 
regimes where the rumor is bounded to a finite neigh- 
borhood of the initially infected site and where it affects 
a finite fraction of the population are separated by a well- 
defined transition at a finite randomness pc of the under- 
lying small-world network. The results of the finite-size 
scaling analysis support the presence of a critical phe- 
nomenon at Pc- Quite interestingly, the transition is also 
found if, instead of building a frozen small-world net- 
work as above, a dynamic small- world is considered [ pT| . 
In this case, each contact of an I-element is established 
with one of its 2K nearest neighbors with probability 
1— p, and with a randomly chosen site with probability p. 
The transition observed for these dynamic small-worlds 
is of the same type as on static small-world networks, 
but presents quantitative differences. For K = 2, for 
instance, the critical randomness shifts to pc ~ 0.07. 

The geometrical properties of small-world networks 
have been shown to exhibit a cross-over from ordered 
to random behavior at p A^~^, which implies a van- 
ishing critical value of p for asymptotically large systems 
|^|-^. In contrast, we have here presented evidence that a 
dynamical process occurring on such structures exhibits 
critical behavior at a finite value of the randomness. This 
difference suggests that an explanation for the origin of 
the transition in purely geometrical terms is unsuitable. 
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and specific dynamical properties must be taken into ac- in Statistical Physics (Springer, Berlin, 1988). 

count. Very recently, evidence of critical behavior at 
finite small-world randomness has been reported for a 
strictly epidemiological model, whose evolution rules are 
substantially different from those considered here p2[ . 
In particular, the epidemiological model allows for the 
transformation R — > S, giving rise to a closed disease cycle 
(SIRS) through recovery. Moreover, the transformation I 
^ R — !■ S is fully deterministic. The critical phenomenon 
found in that case is a transition to global synchroniza- 
tion of local disease cycles, whereas in our propagation 
process we have a kind of percolation phenomenon |^,|lO| . 
In spite of these basic differences, the dependence of the 
respective order parameters on the randomness p is strik- 
ingly similar, even if compared quantitatively. This sim- 
ilarity calls for further investigation in order to identify 
and characterize the whole class of small-world dynami- 
cal processes with critical behavior at finite randomness, 
and to give an analytical description of such behavior. 

Enlightening discussions with M. Kuperman and G. 
Abramson are gratefully acknowledged. 
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